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Abstract 



2 QECC and Pauli Channels 



The calculating of the coherent information is a fun- 
damental step in obtaining the quantum capacity of a 
quantum channel. We introduce orthogonal and com- 
plete code basis to evaluate the coherent information 
per channel use when the input is the maximal mixture 
of stabilizer codewords. In the code basis, the output 
density matrix is diagonal, the joint output of the sys- 
tem and the auxiliary is block diagonal. The coherent 
information is worked out by counting the weights of 
error operators. 

PACS: 03.67.Hk, 03.67.Pp, 03.65.Ud 



1 Introduction 



The basic issue in quantum information theory is quan- 
tum coding theorem. After ten year's efforts, quantum 
coding theorem had at last been proven. The rate of 
faithfully transmitting quantum information per use 
of quantum channel is limited by quantum capacity, 
the capacity is asymptotically achievable [T] [2] [3] [3]. 
Quantum capacity is the maximization of coherent in- 
formation [5] over all input states. Unfortunately, since 
coherent information is non- additive [BJ, quantum ca- 
pacity in single letter form is not available except for 
degradable 7 or anti-degradable channels. Regulation 
is need, that is, block input with infinitive number of 
qubits should be used to calculate quantum capacity in 
general. We may only obtain the lower bound of quan- 
tum capacity. The obstacle of obtaining the coherent 
information with other multipartite input state is obvi- 
ous, the dimension of the state increases exponentially 
with the the number of the input qubits, making the 
calculation of the output entropy and the entropy ex- 
change (thus the coherent information) an awful work. 
We will greatly reduce the complexity of diagonaliz- 
ing the output density matrix by introducing quantum 
error-correcting code (QECC) as the input state. 



The theory of QECCs was established more than a 
decade ago as the tool for fighting decoherence in quan- 
tum computers and quantum communication systems 
[5]. Maybe the most impressive development in quan- 
tum error-correction theory is the use of the stabilizer 
formalism^ [TD] P] [E]. The power of the stabi- 
lizer formalism comes from the clever use of group the- 
ory. The n-fold Pauli operators {/, X, Y, Z}® n together 
with the possible overall factors ±l,±i form a group 
G n under multiplication, the n-fold Pauli group. Sup- 
pose S is an abelian subgroup of G n . Stabilizer coding 
space T is the simultaneous +1 eigenspace of all ele- 
ments of S, T = : M = \ip) ,VM 6 S}. For an 
[[n, k, d]] stabilizer code, which encodes k logical qubits 
into n physical qubits, T has dimension 2 k and S has 
2 n ~ k elements. The generators of S are denoted as Mj 
(i = 1, . . . ,n— k) which are Hermitian. There are many 
elements in G n that commute with every elements of 
S but not actually in S. The set of elements in G n that 
commute with all of S is defined as the centralizer C(S) 
of S in G n . Clearly S C C(S). 

Denote il — Yli(I + M)- Due to the properties that 
the elements of stabilizer group S commute and Mf = 
I, it follows Mj 11^1 + Mi) = Yliil + Mi). Further, we 
have flMjfl — Q 2 = 2 n ~ k fl. For error operator E a that 
anti-commutes with at least one of the generators Mj, 
we have 

nE a n = o, (i) 

this is due to (I + Mi) (I - M t ) = I - M? = 0, the 
operator factor (7 — Mj) comes from E a (I + Mi) = 
(I — Mi)E a if E a anti-commute with Mj. 

In Krauss representation, Pauli channel map E act- 
ing on qubit state p can be written as £ (p) = fp + 
p x XpX + PyYpY + p z ZpZ, where p x (y,z) S [0,1] are 
the probabilities, / = 1 — p x — p y — p z G [0, 1] is 
the fidelity of the channel. For depolarizing channel, 
Px = Py — Pz = P, f = 1 — 3f>. The total error prob- 
ability is 3p. For n use of depolarizing channels with 
n qubits input state p, we have the output state p' = 
£®n {p) = j- a Va E aP El with Va = /" ' * 'p,p:,p[ for 
E a = X'Y^Z 1 . The purification of p is , we have the 
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joint output state p e = (£® n <g> I® n ) (|#) (*|) , whose 
entropy is the entropy exchange. The coherent infor- 
mation is I c — S(p') — S(p e ), where S(-) is the von 
Neumann entropy. 



3 The coherent information of 
depolarizing channel with 
coded input 

3.1 The [[5,1,3]] code 

To illustrate non-zero capacity even for zero fidelity 
of Pauli channel, we first consider the example of 
[[5,1,3]] QECC as the input state to the depolar- 
izing channel. The stabilizer code has four gener- 
ators Mi = X 1 Z 2 Z 3 X 4: , M 2 = X 2 Z 3 Z 4 X 5 , M 3 = 
X 1 X 3 Z 4 Z 5 , Mi = ZiX 2 X 4 _Z 5 . The codewords_are 
|0) = ^ 100000} , |T) = \xn |00000) , where X = 
XXXXX,Ji C C(S)\S. Another useful operator _m 
C(S)\S is Z = ZZZZZ, which anti-commutes with X. 
The input state p is chosen to be ^(|0) (0| + |l) 
with maximal input entropy in the logical qubit ba- 
sis |0) and |l). To evaluate the eigenvalues of output 
state and joint output state of the system and the aux- 
iliary, we introduce the code basis: J) = Aj |0) ( J = 
0_^ . . , 31), where Aj_= /, X m (to = 1, . . . , 5), Y m , Z m 
,X, X m X ,Y m X , Z m X. . Note that [[5, 1, 3]] can correct 
any single qubit error, the basis are orthonormal. For 
input code state |0) , the elements of output density 
matrix are ( J\ £ a VaE a |(J) (0| E\ \K) . Denote |0) = 
100000) , the factor (J\ E a |0) = i (0| ClA 1 jE a Cl |0) is 
when Aj£ a ^ S + SZ (up to ±l,±i factors, which have 
no effect in channel mapping and are omitted hereafter) 
, and is 1 when A^jE a E S + SZ , according to the 
following reasons. We have three cases: (i) A^jE a 

r anti-commutes with some genera- 



diagonalized. 



C(S), then A]jE, 



tors of S, by Eq. (TJJ), so (J\E a |0) = 0; (ii) A],E a E 
C(S)/(S+SZ), equivalents k},E a E X(S + SZ). Then 
there is some M G S such that (0| flA]jE a fl |0) = 

J6 (o| nxn |o> - (o| T) = o ; (hi) 

SZ, then there is some M E S such that 



±(o\nMnx\o) = 



A>jE a E_S 

(J\E a \0) = ^(0|OMfi|0) = i(0|OO|0) = 1. If 
A}jE a G S + SZ and E\A K G S + SZ, we have their 
product A\E a ElA K = A\A K G S + SZ , which is only 
possible when J = K, since our Aj is so chosen that 
each of which is the head of one of cosets G5/ (S + SZ) . 
Thus the state £ ® 5 (|0) (0|) is diagonal in the code ba- 
sis. Similar result can be found for £® 5 (|l) (l|). The 
output state in the representation of the code basis is 



p'jk 



\5jK 



(E 



Va 



(2) 



with the conditions of AjE, 



E S + SZ for the first 
SZ for the second term at the 



term, and XA*jE' a E S 
right hand side . 

The purification of the input state p could be l^) = 
4=(|0~)|0a/ + |I)|1^i/) ; wnere the first logical qubit is 
for the system , the second logical qubit is for the aux- 
iliary and denoted by the subscript a. The joint output 
state of the system and the auxiliary is p e — (£ ® 5 (g> 
I 05 ) (|*) (\fr|). An obvious basis for the joint out- 
put density matrix is | J) 0^4) , | J) lyt) . The non-zero 
elements are (J| (0 A \ p e \ J) \0a) , (J\ (0 A \ Pe \J') |Ta) , 
(J' I (T A | p e I J) |0i) , (J'\ (Ta\ Pe I J') \Ta) , where J' = 
mod(J + 16,32). The reason can be seen from the 
calculating (J\ (0a\ p e \K) \1a) for example. We 
have (J\(0 A \p e \K)_\T A ) = \ (J| ^ 5 (|tT) (T|) \K) = 
I So Va (^1 0) (l| El \K) , which is nonzero when 
G S + SZ and E\AkX G S + SZ, thus 



A]jA K X 



E S + SZ, hence K 



PeJ 



A\E a ElA K X 

mod(J+ 16, 32) = J'. In the basis of | J) \0 a ) , \J) \l a ) 
by rearranging the subscripts, the matrix p e can be de- 
composed to the direct summation of 32 submatrices. 

Pe = ffij=oPeJ, with 

<J| (0a| p e I J) |CU) (J\ (0a\ Pe \J') \Ta) 
(J'\ (1a\ pe \ J) |CU) (J'\(l A \pe\J')\lA) 

We have (J| (0 A Jp e JJ) \0 A ) = | <J|£® 5 (|0) <0|) |J) = 

\Y,aVa{J\E a \Q)(Q\El\J) = ^ E _ aeAj( s + SZ)Va, 

and_(J'|(U| / o_ e |J')|U) = (J\(0A_\pe_\J)\0A); 

(J\(0 A \pe\J')\lA) = kEaVa(J[Ea\0)(l\Et\J'}. 

Notice that (T| E\ \ J') = (0| XE\Aj, |0) = 
(0| XE^AjX |0) , if A\E a E S, we have X~E\Aj~X = 
E\Aj, while for A]jE a G SZ, we have ~XE\AjX = 
-E\Aj. Thus 



PeJ - 7: 



VeJ + VoJ 
VeJ ~ VoJ 



VeJ ~ VoJ 
VeJ + VoJ 



with 



VeJ 
VoJ 



zZ ^ 

E a eAjS 
E a eA.jSZ 



(4) 

(5) 
(6) 



Thus the eigenvalues of p e j are f] e j and rj j. We can 
obtain from @ that the eigenvalues of p' are 

^J=2 ^oJ + VeJ + VoJ' + VeJ')- (7) 
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The average coherent information per channel use thus 
is 



c° Q = c§ 4 = cg 4 = (15,0,42,0,7,0,0,0); (ii) 7 fold 
degenerate represented by cf = (6,24,21,8,4,0,1,0), 



_ ! 32 

Ic = - ^(-Aj log 2 Aj + r) j log 2 r] j + 7y eJ log 2 r/ eJ ). 
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(11,16,18,16,3,0,0,0), c; 
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J=0 



(8,19,24, 10, 0,3,0, 0);(iii) 7 fold degenerate rep- 
resented by C | = cf, eg = c§ 5 , c? 2 = cf 2 = c?;(iv) 



(8) 7 fold degenerate represented by cf 5 = cf, C 15 
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The eigenvalues r? e j ,?7 j can be obtained by count- 
ing the weights of all the operators of AjS and AjSZ, 
respectively. The weight of an error operator is the 
number of qubits on which it differs from the iden- 
tity. An error E a G AjS with weight j will contribute 
rj a = f n ~ip> to the eigenvalue of r] e j, similarly, An 
error E a G AjSZ with weight j will contribute r] a = 
f n ~ J P J t° the eigenvalue of r) j. Thus r\ e j and ?7 jcan 
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c§ 5 ; (v) 42 fold degenerate repre- 
sented by c| 3 = c° 3 = c% 7 = (8,21,20,10,4,1,0,0), 
c° 7 = (9,20,18,12,5,0,0,0). The eigenvalues of the 
output state p' are Aj = \{r) oJ + rj eJ + r] j> + i] e j>), 
with J' = mod(J + 64,128). The average coherent 
information per channel use then is 



be written as V; cf /" •</-•< and V; Q cff n 'p> T c = i[- 
,respectively. Furthermore, they are characterized 
by vector Cj = (c§ J , c| J , dj J , c| J , J , cg J ) and Cj = 
(cg J ,c| J ,cg^,c^ J ,cf 7 ,cg J ). A detail counting shows 
eg = (0, 15, 0, 0, 0, 1), eg = cf 6 = cf 6 = (6, 0, 10, 0, 0, 0). 
The eigenvalues of the subscripts J = 1,2,3,4,5, 
22, 23, 24, 25, 26, 11, 12, 13, 14, 15 are degenerated and 
represented by cf = (3,8,4,0,1,0) and cf = 
(4,6,4,2,0,0); also the eigenvalues are degenerate for 
J = 17, 18, 19, 20, 21, 6, 7, 8, 9, 10, 27, 28, 29, 30, 31 and 
represented by cf 7 = cf 7 = cf . The average coherent 
information per channel use then is 



1 r (jho + 3t? o0 ) (?7 e o + 3t7 o0 ) 
j[ 2 °S2 2 

+VeO log 2 r/ e0 + 3f7o log 2 r] o0 ] 



(r] el + 2t] o1 + l] e65 ) (?? e i + 2?7 i + 77 e65 ) 

+6[ 2 g2 2 

+Vel log 2 J? e l + 2?7ol log 2 T] ol + f] e65 log 2 77 e65 ] 
LR[ (^7e87 + 3?7o87) (>7e87 + 3?7 87) 

~ 2 log 2 - 

+Ve87 l0g 2 Ve87 + 3r? o8 7 log 2 ^087] ■ (10) 



Ir = 



1 r {VeO + 3t7 o0 ) 1 (?7eQ + 3?? o) 

5 [ 2 g2 2 

+7y e0 l0g 2 T] e0 + %VoO log 2 7? o] 
, q r (??el + 37? i) (Tfei +3?7 Q i) 

+3[ ^ log 2 - 

+»7el l0g 2 %1 + 3r/ i log 2 T] ol ] . 



(9) 



3.2 The [[7,1,3]] code 

The six stabilizer generators of [[7,1,3]] code 
are M l = X X X 2 X Z X^ M 2 = X x X 2 X h X % , 
M 3 = X^sXsXj, M 4 = Z 1 Z 2 Z 3 Z 4 , M 5 = Z 1 Z 2 Z 5 Z 6 , 
Mq = Z1Z3Z5ZJ. The bit flip and phase flip operators 
for the encoded (logical) qubit arc X = X 5 X e Xr 
and Z = Z^ZqZj, respectively. Since the code 
can correct any single qubit errors of X, Y, Z 
types and meanwhile it can correct errors XiZj 
{i j) type, it is convenient to choose the 

128 coset heads as Aj (J = 0, . . . , 127) = I, 

X k ,(k = l,...,7),Y k ,Z k , XiZj i,j = 1 7), 

X, X k X,(k = 1,...,7), YfcA, Z k X, XiZjX (1 ^ 3: 
i,j — 1,...,7). With i] e j and rj j being written 

as V' „«•'/'/" J l> ! and Ej=o c f/""V' the non- 
zero eigenvalues of the joint output state p e now 
can be characterized by vectors Cj and Cj. They 
are (i) non-degenerate eg = (0,42,0,21,0,0,0,1), 



3.3 The [[8,3,3]] code 

The improvement to the lower bound of quantum ca- 
pacity comes from [[8,3,3]] code. The code encodes 
3 logical qubits into 8 physical qubits and corrects 1 
error. Thus the coding rate is higher than [[5,1,3]] 
code for noiseless channel. And it indeed provides 
a tighter lower bound for the quantum capacity at 
high fidelity domain. [[8,3,3]] code has generators 
Mi (i = 1,...,5) of its stabilizer group S , with 

Mi - nti *<> M 2 = nti Z U M 3 = X 2 X 4 Y 5 Z e Y 7 Z s , 
M A = X 2 Z 3 Y i X 6 Z 7 Y 8 , M 5 = Y 2 X 3 Z 4 X 5 Z 6 Y 8 . In 
the centralizer C(S) there are bit flip and phase flip 
operators Xj, Zj (j = 1,2,3) for the encoded (logi- 
cal) qubit, with X~i_= A"iA" 2 Z 6 Z 8 , Xj. = X^Z^Zj, 
X3 = X1Z4X5ZQ, Z\ = Z 2 Z4ZgZg,, Z 2 = Z^ZjZ jZ^, 
Z% = Z^ZqZ'jZs,. The codewords are \k\k 2 kz) = 

j^X^X^X^fl \0) , with ki = 0, 1. The channel in- 
put state p could be chosen as the equal probabil- 
ity mixture of codeword states, each codeword has 
a probability of |. The basic set of the heads of 
cosets G$/(S x Z) (Z is the group with generators 
Zj) can be chosen as the correctable single qubit er- 
rors X m , Y m , Z m and some other two qubit errors such 
as XiX i+ i (i = 1, . . . , 7) and the identity /. The num- 
ber of the elements in the basic set is 32. The whole 
coset head set is obtained by the multiplication (at 
right) of the basic set with group X whose genera- 
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tors are Xj. Denote the elements of the coset head 
as Aj (J = 0, ...,255). Suppose the input state is 
P = iE* 1 ,fe,fc 3 =o | ^1^2^) (fcifofoli then the purifica- 
tion of pis i*) = ^e!,,),^,^ |^ifc 2 fc 3 ) |(fcifc 2 fc 3 )A) , 

the joint output of the system and auxiliary is p e = 
(£® 8 ® I® 8 ) (|*) (#|) . In' the basis of \J) = Aj |000) , 
we have 



& = ELoCij+i/ n "V with 



c = 



(^|((W^ A |Pe|J)|(ili2i3)A) 

^ % (A-| £ a |fcK) (ju^| St | j) , (ii) 



which is nonzero when Iiljl* 3 e 
5 x Z and ^AjXf X^xf e S x Z, 
thus X^X^^A^K^A^X?^ 3 

X^X^X^A^AjXf Xjxf e 5 x Z. How- 
ever, Aif and Aj are coset heads, so we 
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have A 



K 



-r^k 1 +j 1 —k 2 +j2^rTk 3 +j. 



L x 



X; 



Let 



X^ x\ 2 x\ 3 K^ K E a = SbZl 1 Z c 2 2 Zf , where S b 
(6 = 0,..., 31) is the element of stabilizer group 5* and 
d = 0,1. Then E a = AkX^X^X^ 1 S^Z^T 3 = 

We 



The degeneracy vector is d =(e?i, . . . , di4)=(l, 35, 
28,112, 168,56, 112,24, 504,168, 168,336, 168,168). 
The entropy exchange is 

14 

S (Pe) = -J2 d ^ l °S2^i)- (13) 

The eigenvalues of the output state p' are 
1 



have 



(000| tf K \ K X**X*fxfs h T*T*^X^X*fX^ 1 000) = 

(_ 1 )fe 1 c 1 + fe 2C2 +fe 3 C3 and (jlj 2 j3\El\J) 

(000| X^X^Xf ZfZzZl 1 S ] b X^X^X j3 A^Aj |000) = 
(_ 1 )jici+j 2 c 2 +j 3 c3_ So that 



Ai 
A 2 
A3 

A 4 

A 5 



g tfi +356 + 28^), 
4(6 + 6), 

86 + 76), 



^4 



E?=i(fe.+J.)c 



with 



i(6 + 6 + i46), 

^[6 + 7(36 + 6o 
+61 + 262 + 6s + 64)] 

vector d' 



8d', 



(14) 
and 



Va>- 



(12) 



degeneracy < 

d'=(d' 1 ,...,4)=(l,3,2,2,24). Thus the entropy 
of p' is 



Where E a , e A 7 Xf S 
and Tnod(K + 32(4fc 3 + 2fc 2 + fci),256) = 
mod(J + 32(4j 3 + 2j 2 + ji),256). By rearranging 
the basis, the joint output state p e can be writ- 
ten in a block diagonalized form with each block 
being a 8 x 8 submatrix. A detail analysis shows 
that each 8x8 submatrix can be diagonalized 
with Hadamard transformation. The eigenvalues 

of p e are ^aeAjS^a, EaGAjSZi Va, J2a£AjSZ 2 Va, 
EaeAjSZiZa Va, Sa6A. 7 SZ 3 Sa£AjSZiZ 3 ^O' 

EaeA jS z 2 z 3 Va, EaeA., sz{z 2 z 3 Va- The total number 
of the nonzero eigenvalues of p e is 256 x 8 = 2048. 
A detail counting shows that each of the eigen- 
value should be on of the & (i — 1,...,14), where 



S(p') 



i^^log^A;). 



(15) 



»=i 



The average coherent information per channel use is 

5 1 14 

T ° = ~ E + o E d ^ lo S2te)- (16) 



i=l 



? = 1 



The coherent information is shown in Figure 1 and Fig- 
ure 2. 



4 General Pauli channels 

For Pauli channel with 3 error probabilities p x , p y , p z , 
the coherent information of coded input can be eval- 
uated in the same way as depolarizing channel. The 
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Figure 1: Solid line for [[5,1,3]]; Dash line for [[7,1,3]]; 
Dot-dashed line for [[8,3,3]]. 
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Figure 2: Solid line for [[5,1,3]]; Dash line for [[7,1,3]]; 
Dot-dashed line for [[8,3,3]]; Dot line for hashing 
bound. 



Table 1: The eigenvalues of p e 
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only difference is that now we should count the num- 
bers of the each type of error separately in calculating 
^2 a r] a . Define the functions 

gi(t,z,x,y) = t 5 + 5t(x 2 y 2 +x 2 z 2 + z 2 y 2 ), 
g 2 {t,z,x,y) = t 4 x + 2t 2 x(z 2 + y 2 ) 

+2tzy{2x 2 + z 2 + y 2 ) 
+x(x 2 y 2 + x 2 z 2 + z 2 y 2 ). 

Then the eigenvalues of the joint output state p e are 
given in Table 1, 

where the eigenvalues are expressed as with de- 
generacy di. The entropy exchange then is 

4 4 

S iPe) =~J2 di J2 & l0 §2 & ■ ( 17 ) 
i=l j=l 

The eigenvalues of p' are Xi = \ Y^=i £y with degen- 
eracy 2di, thus we have 

4 

^(-°') = -^2d l A 4 log 2 A l . (18) 

i=l 

The coherent information per channel use for Pauli 
channel with [[5, 1,3]] code as input is 

4 4 4 

t c = -^E^I^^E^) 

t=l j=l 3=1 
4 

3 = 1 

5 Dicussions and Conclusions 



There are bounds on the coding rate of QECC, the 
quantum Hamming bound 13 , Knill-Laflamme (quan- 
tum Singleton) bound |14] , Gottesman bound and so 
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on |10) . The first two give rather tight upper bounds on 
some of additive quantum codes. The quantum Ham- 
ming bound (hashing bound) is a strict upper bound 
for non-degenerate (pure) quantum code, as it is seen 
from figure 2, where all three average coherent infor- 
mation calculated are upper bounded by the hashing 
bound. However, it has been known that quantum 
Hamming bound can be violated by degenerate (im- 
pure) quantum codes [BJ. They obtain the result by 
calculating the coherent information of depolarizing 
channel with repetition quantum codes. We have in- 
troduced a systematical way of calculating the coher- 
ent information of Pauli channel with quantum code 
as input state. The main finding of this paper is that 
the channel output density matrix as well as the den- 
sity matrix of the joint output of the system and the [14] E. Knill, R. Laflammc, Phys. Rev. A 55, 900 
auxiliary can be diagonalized for Pauli environment (1997). 
with quantum code as input, the eigenvalue problem 
is reduced to counting the weight of the error oper- 
ators in the coset. We have presented the input of 
[[8,3,3]] code as an example of calculating the coherent 
information of input state with multiple logical qubits. 
It is anticipated that our method should promote the 
way of violating quantum Hamming bound by calcu- 
lating the coherent information of depolarizing channel 
with quantum code of encoding multiple logical qubit. 
Meanwhile, our method provide the way of calculat- 
ing the lower bound for the distillable entanglement 
of quantum code state passing through Pauli channel, 
according to hashing inequality [3] ■ 

Funding by the National Natural Science Foundation 
of China (Grant No. 60972071), Zhejiang Province Sci- 
ence and Technology Project (Grant No. 2009C31060) 
are gratefully acknowledged. 

References 

[1] S. Lloyd, Phys. Rev. A 55, 1613 (1997). 

[2] H. Barnum, M. Knill and M. A. Nielsen, IEEE 
Trans. Inf. Theory, 46, 1317 (2000). 

[3] M. Horodecki, P. Horodecki and R. Horodecki, 
Phys. Rev. Lett. 85, 433 (2000). 

[4] I. Devetak, IEEE Trans, inf. Theory 51, 44(2005). 

[5] B. Schumacher and M. A. Nielsen, Phys. Rev. A 
54, 2629 (1996). 

[6] D. P. DiVincenzo, P. W. Shor and J. A. Smolin, 
Phys. Rev. A 57, 830 (1998). 

[7] I. Devetak, and P. W. Shor, Comm. Math. Phys. 
256, 287(2005). 



[8] P. Shor, Phys. Rev. A 52, R2493 (1995). 

[9] A. R. Calderbank, E.M. Rains, PW. Shor, N. J. 
A. Sloane, Phys. Rev. Lett. 78, 405 (1997). 

[10] D. Gottesman, thesis, California Institute of Tech- 
nology, Pasadena, CA (1997). 

[11] D. Gottesman, Phys. Rev. A 57, 127 (1998). 

[12] M. A. Nielsen, I. L. Chuang, Quantum Computa- 
tion and Quantum Information (Cambridge Univ. 
Press, New York, 2000). 

[13] A. Ekert and C. Macchiavello, Phys. Rev. Lett. 77, 
2585 (1996). 



6 



